FERMIONIC ENTROPY of the VORTEX STATE in d-WAVE SUPERCONDUCTORS. 
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Abstract 



In the d-wave superconductors the electronic entropy associated 
^J ■ with an isolated vortex diverges logarithmically with the size of the 

~f*. . system even at low temperature. In the vortex array the entropy per 

0> vortex per layer, Sy, is essentially larger than k,B and depends on the 

distribution of the velocity field v s around the vortex. If there is a 
£— j . first order transition with the change of the velocity distribution, then 

r^" there will be a big entropy jump ASy ~ Iks at the transition. This 

entropy jump comes from the electronic degrees of freedom on the vor- 
tex background, which is modified by the vortex transition. This can 
explain the big jump in the entropy observed in the so-called vortex- 
melting transition |lj, i n which the vortex array and thus the velocity 
field are redistributed. The possibility of the Berezinskii-Kosterlitz- 
Thouless transition in the 3- dimensional d-wave superconductor due 
to the fermionic bound states in the vortex background is discussed. 
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1 Gap nodes and scaling. 

The low energy properties of the superconductors with nodes in the energy 
gap are governed by the electronic excitations close to the gap nodes. The 
electronic density of states (DOS) in the homogeneous superconductor is (see 
eg the Review paper 0) 

N(E) ~ N F (Jp) , £«T C , (1) 

where D is the dimension of nodes and Np is the DOS on the Fermi level 
in the normal metal. In the mixed state of superconductor the superflow 
around the vortex leads to the Doppler shift of the energy E to E + v s (r) • p, 
where v s (r) is the local superfluid velocity. This gives the finite DOS at zero 
energy 

N(0)~N F (^y~ D . (2) 

Here v s is the characteristic value of the superfluid velocity in the vortex 

/o N.l/2 

array: v s ~ h/m 3 Rv } where Ry ~ £ ( -^ J is the intervortex distance, B is 
magnetic field and £ ~ Vp/T c is the coherence length. Thus Pfv s /T c r b 



. B ^ 
and this gives the following electronic DOS at zero energy as a function of 

magnetic field 

, B . l-D/2 

N(Q) - N F (-J . (3) 

There are two different regimes of strong and weak fields with the crossover 
parameter 

T T (B c2 \ l l 2 

x ~k^rTA-B-) ■ (4) 

which separates the superflow dominating regime j;<1 from the temper- 
ature dominating regime x ^> 1 (see M). In general the thermodynamic 
functions depend on the parameter x. For example, the free energy of the 
excitations on the background of the vortex array is 

F(T,B)=N F T^—J F(x) (5) 

where F(x) is the dimensionless function of the dimensionless parameter x 
(see also the recent paper by Simon and Lee |5[]; their crossover parameter 



differs from our by the factor JT c jEp. This is because Simon and Lee used 
the linearized spectrum of the fermions in the very vicinity of the gap nodes, 
which can be justified only at rather low temperature, T <C T^jEp). The 
normalization can be found from the low field asymptote x >> 1, where the 
largest contribution comes from the bulk superconductor, while the effect of 
vortices is small. It follows from Eq.(l) that the free energy of the homo- 
geneous state is ~ — N F T^{T /T c2 ) A ~ D , and this gives the normalization of 
F(T, B) and the low field asymptote 



F(x) ~ -x"~ D , x > 1. (6) 



2 Scaling for d-wave superconductor. 

2.1 Free energy 

In the case of nodal lines, ie for a node dimension D — 1, one has the 

following estimate for the free energy of the excitations in the background of 

the vortex array: 

/ B \ 3/2 - 
F(T,B) = N F T? — F(x) (7) 



B, 



r-l 



Let us find the asymptotes of F[x) at x ^> 1 and x<l. Let us consider 
first the weak-field case x ^> 1, i.e., the case of a dilute vortex array. In 
addition to the largest asymptote F(x) ~ — x 3 from the bulk superconductor 
there is also the contribution from vortices. It comes from the modification of 
the normal- component density due to excitations: p n {T) ~ pN(T)/N F . This 
leads to a decrease of the kinetic energy of superflow around the vortex, thus 
the contribution to the energy of the vortex array with the vortex density 
n = B/<&o, where $o is the flux quantum, is 

F(B - 0) = -\np n {T) f d\ v s 2 ~ -N F T?JL?- In ^ . (8) 

2 JR v >r>hv F /T B C 2± C £-tc 

This corresponds to the term — a; In a; in F(x). So, the leading terms in the 
asymptote of F(x) at x ^> 1 are 

F{x) ~ — xhix — x 3 , x ^$> 1. (9) 



The asymptote of F(x) at x -C 1 is 

F(x) ~ -1 -x 2 , x< 1 (10) 

Here both terms are from the vortices. The first one is temperature inde- 
pendent and does not contribute to the entropy or specific heat, but does 
contribute to the vortex magnetization. It comes from the nonzero normal- 
component density at T = due to the superfluid velocity: p n {T = 0) ~ 
P N(E = 0)/N F ~ p( PF v s /T c ) 2 ~ D (see § g g 1): 

The second (quadratic in x) term in Eq.(10) gives a term linear in tem- 
perature to the specific heat: C(T, B) oc T\/B~ ||. 

2.2 Vortex entropy 

For the entropy density one has 

^,5) = -^ = ^— T c 5(x) (12) 

OF 
S( X ) = -_ (13) 

The asymptotes of S(x) at x ^> 1 and i-C 1 are 

5(ar) ~ In x + x 2 , x > 1 (14) 

5(x) ~ x , x < 1 (15) 

The vortex part of the normalized entropy S(x), i.e., without the bulk term 
x 2 in Eq.(14), can be written using the interpolating formula 

S(x) ~ ln(x + 1) (16) 

which gives both the logarithmic term in Eq.(14) at large x and linear term 
in Eq.(15) at low x. 



It is instructive to write the vortex entropy per vortex per layer: 

S V (T,B) E F 

— ln(ar + l), (17) 

k B T c 

where Ep is the Fermi energy. Note that the logarithmic vortex entropy 
in Eq.(14) also follows from the 1/E behavior of the vortex DOS found in 
Ref. ||]. Using the result of Ref.Q one can find an exact equation for the 
vortex entropy at large x using an axial distribution of the superfluid velocity 
around the vortex, v s = [Tijlm^r)^: 

^^ = 21n2^1nz , x » 1. (18) 

k B A' v ; 

Here A' is the angle derivative of the gap at the node. 

2.3 Heat capacity. 

For the heat capacity one has 

C(T,B) = T— = N F T C —C(x) (19) 

C{x) = x- (20) 

Using the interpolating formula (17) for the entropy, one obtains interpolat- 
ing formulas for the vortex part of the heat capacity: 

C(T,B)~N F T*-?-, (21) 

-D C 2 1 + X 

which gives both asymptotes found in Ref. 0: C(x) ~ s for i C 1 and 
C(x) ~ 1 for x > 1. 



3 Discussion. 

The electronic entropy per vortex per layer Sy in Eq.(17) is much larger 
than ks even at T <C T c . For an isolated vortex this entropy diverges as the 



logarithm of the dimension R of the system: Sy ~ kB{EF/T c ) \nR (actually 
R is limited by the penetration length) and is at least a factor of EpjT c ^> 1 
larger than the configurational entropy of the vortex in the 2-dimensional 
system, S^onf ~ fc^hi-R. The logarithmic behavior of Sy, with R limited by 



the intervortex distance Ry, persists till T ~ T C JB/B C 2 (or x ~ 1). Due 
to the large factor Ep/T c ^> 1 the entropy per vortex per layer can be of 



order kg even at T < T C J B / ' B c2 (or x < 1), but it finally disappears in the 



high-field limit (T < T C ^B/B c2 or x < 1). 

It is important that Sy depends on the distribution of the velocity field v s 
around the vortex. If there is a first-order transition with the change of the 
velocity distribution, one can expect a big entropy jump, AS ~ Iks- This 
entropy jump comes from the electronic degrees of freedom in the vortex back- 
ground, which is modified by the vortex transition. This can explain the la- 
tent heat L ~ 0.45fc#2"y vortex/layer ovserved on the so called vortex-melting 
line in a detwinned Y-123 crystal |1(| and L ~ 0.6 ± Cl^eT/vortex/layer in 
a twinned sample of Y-123 0]. Even higher values of the entropy jump have 



been deduced from the magnetization measurements [pjj| . Such an entropy 
jump can occur both at the vortex-melting transition and at a first-order 
transition in which the structure of the vortex lattice changes, say, from a 
hexagonal lattice close to T c to a distorted tetragonal lattice far from T c . The 
latter structure was observed in Ref . [jnj and discussed in Ref . [|T3|] . 



Note that the fermionic entropy of the 3-dimensional vortex loop of the 
length R is oc R In R, as distinct from the configurational entropy of the loop 
S con { oc R. This R In R behavior of the vortex loop entropy, together with 
the large prefactor, can in principle cause the Berezinskii-Kosterlitz-Thouless 
transition in the 3-dimensional system. This is supported by the following 
observation. It appears that if one uses the quantum-mechanical approach to 
the vortex DOS, by calculating the discrete bound states of the fermions on 
the background of the inhomogeneous distribution of the superflow around 
the vortex, one obtains a value of the vortex DOS that is twice as large 
as that obtained from a classical treatment of the fermions in terms of the 
Doppler shifted energy E + v s (r) • p (Ref. [Q). 

The classical approach gives the conventional expression for the energy 



of isolated vortex in terms of the superfluid density p s (T) = p — p n (T) 

Fv = ^(P - Pn(T)) In | . (22) 

At low T, where p n (T)/p ~ T/T c , the second term corresponds to a loga- 
rithmic entropy of the vortex. The quantum-mechanical approach in terms 
of the fermionic bound states in the vortex background suggests the larger 
contribution of the fermions to the vortex entropy. This can be written using 
the enhancement factor 1 + a(T) 

F v = ^[p-(l + a(T))p n (T)]lnj . (23) 

According to [§J one has a(0) = 1. If this value of a persists to higher 
temperatures, the energy of an isolated vortex becomes zero at some tem- 
perature Ty < T c , where p n (Ty) = \p. Thus at Ty one would have the 
Berezinskii-Kosterlitz-Thouless transition in the 3- dimensional system, oc- 
curing due to the essential contribution of the fermionic bound states to the 
vortex entropy. However, it is more natural to expect that a(T) decreases 
continuously with temperature, approaching zero value at T c , since the effect 
of the bound states should be negligible in the Ginzburg-Landau region. So, 
the possibility of the Berezinskii-Kosterlitz- Thouless transition in this 3D 
system depends on the details of the behavior of a(T). 
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